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Ultrashort non-resonant electromagnetic pulses applied to effective one-electron systems may operate on the 
electronic state as a position or momentum translation operator. As derived here, extension to many-body 
correlated systems exposes qualitatively new aspects. For instance, to the lowest order in the electric field 
intensity the action of the pulse is expressible in terms of the two-body reduced density matrix enabling thus to 
probe various facets of electronic correlations. As an experimental realization we propose a pump-probe scheme 
in which after a weak, swift "kick” by the non-resonant pulse the survival probability for remaining in the initial 
state is measured. This probability we correlate to the two-body reduced density matrix. Since the strength of 
electronic correlation is bond-length sensitive, measuring the survival probability may allow for a direct insight 
into the bond-dependent two-body correlation in the ground state. As an illustration, full numerical calculations 
for two molecular systems are provided and different measures of electronic correlations are analyzed. 

PACS numbers: 82.53.Kp,31.10.+z,31.15.A-,71.10.-w 


I. INTRODUCTION 


Recent experiments demonstrated the feasibility of ultra- 
short attosecond laser pulses and their use in combination with 
other optical pump-probe techniques to access various facets 
of the electronic and ionic dynamics HHS. Of a special rele¬ 
vance here are strongly time-asymmetric pulses. The simplest 
example is a mono cycle, time asymmetric pulse for which 
one half cycle could be short and strong and the other half 
cycle could be longer but weaker in a way that the time in¬ 
tegral over the electric field amplitude vanishes, as required 
for a propagating pulse. When interacting with matter, the ef¬ 
fect of two half cycles could be quite different lEHSl. For 
instance, if the first half cycle (with duration r) is shorter 
than the time of relevant transitions taking place in the sys¬ 
tem, whilst the second, weaker half cycle is much longer, then 
the effect of the whole pulse is mainly governed by the first 
half cycle of the pulse (HCP). Henceforth we refer to this sit¬ 
uation as the HCP case. The second half cycle may act as an 
off-set DC tail (if weak/long enough). In this limit it is con¬ 
venient and sometime sufficient to introduce theoretically the 
notion of HCP kick as a useful idealization of the action of the 
whole pulse ll6l [T 0 H 2 TII : Namely, in the HCP case, the system 
is insensitive to the details of the pulse temporal shape. This 
allows relating the action of HCP to that of a kick, i.e., a 6 
function in time. Substantial simplifications of the triggered 
quantum dynamics follow then. We will demonstrate below 
that by performing measurements in an appropriate setup on 
the system immediately after such an excitation allows to in¬ 
fer the strength of the electronic correlation in the system and 
its dependency on the internal structure, and to quantify the 
entanglement. The experimental arrangement that we suggest 
is the following: We start from the ground state, of say some 
molecular structure (another stationary case is also possible) 
and ’’shake” the electronic system by a weak but short HCP. 
As the system evolves and possibly dissociates we measure. 
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e.g., via the time-resolved photoemission fTH the initial-state 
remaining, i.e. what we will call the survival probability. De¬ 
tecting only the kinetic energies of the dissociating ions upon 
the kick, one may access information on the ionic time evolu¬ 
tion that can be used when measuring the time-resolved sur¬ 
vival probability as a function of the molecular bond (since at 
this particular time one has information on the molecular bond 
distance). Currently we are not aware of any corresponding 
experiment to this proposal, however, the ingredients of the 
suggested setup, such as HCP, time-resolved photoemission 
were demonstrated Il23l425]l . From the survival probability 
we deduce features akin to the reduced two-particle density 
matrix. 


Specifically, the envisaged experiment should measure 
upon the kick the initial state occupation probability as a func¬ 
tion of the evolving internuclear distance. It is this quantity 
which we calculated numerically and will be discussing be¬ 
low. Knowing (from theory or experiment Il26l|22l molecular 
bond distance upon the non-resonant pulse excitation we can 
so image the time evolution of the survival probability. We re¬ 
late it below to the reduced two-particle density matrix which 
vanishes for the single determinantal states and, thus, can be 
used as a measure of electronic correlations. 


So the key point of this study is how to measure and quan¬ 
tify electronic correlations. In Sec. we give an explicit 
expression for the probability of the system to remain in its 
ground state after the application of a d-like pulse. While 
for one-electron systems the quantum dynamics was discussed 
by many authors the nontrivial part tackled here is the many- 
body nature of the problem. In Sec. Ill we put the discussed 
probability in the context of other proposed measures of elec¬ 
tronic correlations: the Frobenius norm of the second cumu- 
lant matrix and the von Neumann entropy. Stretching molec¬ 
ular bonds incur in general a change in the amount of correla¬ 
tions in the electronic subsystem. The proposed experiment is 
illustrated by numerical simulations which show what a pos¬ 
sible outcome may be expected (Sec.|IV[). The dependence of 
different correlation measures on the interatomic distances is 
studied, and analytical expressions in the asymptotic regime 
are obtained. Additionally, our calculations allow to test nu- 
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merically for important inequalities for the entropy measures. 


II. THEORY 

Consider an A^-electron system in its ground state. Its 
electronic properties are completely described by the many- 
body wave function i/ro(l,2,... ,N). Although the choice of 
the gauge is irrelevant for the present discussion we will as¬ 
sume here the light-matter interaction to be given in the length 
gauge by the dipole operators where a denotes a projection 
(determined by the pulse polarization). The coupling opera¬ 
tors are of one-particle type, Hermitian and of 

vector character. By choosing a suitable one-particle basis we 
can cast the light-matter interaction in the second-quantized 
form 6H(t) - - Y: ,• where £“ are the com- 

ponents of the electric field vector. Acting on the system with 
a pulse that has the time-dependence E“if) - E^S^it) with a 
peak amplitude at f = 0 , we can write the state of the system 
at f = e shortly after the pulse as 

= e' i/ro = lAo ■ ( 1 ) 

df is a regularized version of the mathematical d-function on 
[ 0 , e] interval and Sij - - fj dlYadj^jEgd^it) has a dimen¬ 
sion of the action (we use atomic units throughout the text). 
Mapping the pulse generated state on the ground (initial) state 
ifro, i-e., taking the overlap (i/^o|iA+) of Q we obtain the sur¬ 
vival (recurrence) probability to be in i/^o after the pulse. This 
quantity is central to the following discussion. 

The quantum-mechanical average of the operator S — 
YijSijcJcj can be treated like any statistical average using 
the cumulant expansion: 

i f- P 

{^o\^+) = exp( —Si -H —S 2 H- —S 3 -H ■ ■■}. (2) 

The correlation functions are well known; 

Si=<5), (3a) 

S 3 = <5^)-3<52)<5)h- 2<5)1 (3b) 

These averages can be computed in terms of reduced density 
matrices (RDM) 


'D'. = <iAo|cJc;|i/ro>, (4a) 

^Df, = {,^o\c]clcicj\i/co), (4b) 

= {l^/o\c]clcicnClCj\lJ/o}. (4c) 

In order to do so we need some additional notations for oper¬ 
ators. 5" is a «-body operator, i.e. it is given by an expression 
containing n creation and n annihilation operators. Let [5”] 
denote a one-particle operator which is given by the nth power 
of S operator in the first quantization, i.e. 


In these notations: 


Q2^-{S ®S) + {[P])-{[S]f, 

S 3 = -{S ® 5 )- 3 <[ 52 ]g, 5 ) + 


where the M ® ‘^B denotes the normal form of the product M 
and the p- and ^-body operators, respectively. It is defined 
as follows: 


'J 

kJ 

'’■4®'b = ZZ4»M,--44--4 


■0 k,l 



In these expressions i, j, and k, 1, are the p and q'-dimensional 
vectors of indices. 

For our discussion it is instructive to introduce the cumulant 
density matrices Il28ll29l ^A which allow to decompose the p- 
RDM in terms of correlated (connected) p-particle correlator 
and products of lower order correlators l(30l : 


-'D' = 'A', 
1! J r 

1 


= Id', a ‘Df A 'D™ + 3 ^A^ A 'D™ + ^A^™, 


(5a) 

(5b) 

(5c) 


where A denotes the wedge product ED (for mathematical 
details we refer to a treatise on differential forms ll^ where 
A appears under the name exterior product). We have, for 
instance; 


'A'. A ‘Af = i (^A''Af - . 

The 1/n! prefactor on the rhs of Eqs. ([^1 appears naturally 
when the density matrices are written in the first quantization: 


—"D(l,...,n; = f i/rUl,... ,n,n- 1 -1,... ,N) 

n\ J 

X ■ ■ ■ ,n',n -t- 1,... ,N)d(n -H 1,... ,N), (6) 

where i = (r,, i,) denotes a collection of space and spin coor¬ 
dinates. 

Returning back to Eq. Q we find by direct comparison 
.^Si = YjijSi/D‘j. It is, however, the second cumulant that 
gives the lowest order (in |Eo|^) contribution to the survival 
probability, namely 

KiAo|iA+)l^ ~ exp[-S 2 ] ss 1 - S 2 , 

S2 = 25,y5«2D''5 + cr2. (7) 

ijkl 


[5"] = 

ij 


Thus, ©2 can be written as the averaged value of a two-body 
operator (first term in Eq. Q), whereas cr^ = ([5^]) - ([5])^ 
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is computed from 1-RDM. This equation can be written in an 
alternative form. Consider the natural orbital basis (i.e. a basis 
in which 1-RDM is diagonal); 


©2 = ^ SijSki pAf + Mkidijdk, - Sudkj)] 


ijkl 


+ ^(SimSmifi - S iiS mm fifm) 


= 2 SijSu^, + Y, \SimfMl - fm), (8) 

i jkl im 


where /) is the occupation number of the ith natural orbital. 
Despite the fact that the second term in Eq. ([8]l has a form of 
the Fermi golden rule ll^ the latter is only valid for suffi¬ 
ciently long pulses (adiabatic switching). 

Eq. Q is a quite remarkable result as it allows to express 
the averaged 2-RDM in terms of experimentally measurable 
quantities; the survival probability and the mean square de¬ 
viation of an excitation operator. Alternatively, one can use 
the form ([^ to access the cumulant as a correction to the 
single-particle result. A possible experiment that we sketched 
in the introduction, could be a pump-probe setup in which the 
system is excited non-resonantly by a half-cycle pump pulse 
and the probe laser pulse is used to monitor the ground state 
occupation. We envisage an application to molecular sys¬ 
tems in which the amount of electronic correlations and, thus, 
([5 ® 5]) are driven by changes of geometry. On the other 
hand, the one-body part cr^ is a quantity that is weakly depen¬ 
dent on the geometric configuration (see also the discussion 
below) and can also be easily computed ED. 

Eor a faithful correlation measure it is desirable to minimize 
the artifacts coming from the dependence of S on the system’s 
geometry. Consider a system subject to stretching. The values 
of the matrix elements S ij can be quite large as compared to 
the equilibrium geometry. Correspondingly, the value of S 2 
obtained according to Eqs. (|7]), (|^ will be at variance with the 
equilibrium value - the effect that is not necessarily reflecting 
properties of ^A. It is possible, however, to suppress to some 
extent the large contributions to S 2 originating from the diag¬ 
onal matrix elements of S. This can be achieved, for example, 
by performing the experiment with oriented systems and ap¬ 
plying the HCP field in the direction for which (S) = 0. Such 
condition is always possible to achieve for systems which con¬ 
tain Cs as their symmetry subgroup and is quite common. Un¬ 
der this condition also cr^ is weakly dependent on the geome¬ 
try. 

Is it possible to devise a measurement that exclusively 
probes the cumulant density matrix ^A? Our answer to this 
question is negative, based on the fact that the generating func¬ 
tional for the cumulant and the reduced density matrices 

is different from the bosonic-like generator of the evolution 
operator Q. In other words, in order to probe the fermionic 
RDMs one needs a direct coupling to fermionic degrees of 
freedom as Grassmann variables 7* in Eq. Q realize. 


Nonetheless, we will demonstrate below using two numer¬ 
ical examples that the survival probability is a versatile mea¬ 
sure of electronic correlations and will compare it to other pro¬ 
posed measures of electronic correlations and entanglement. 


III. MEASURES OF ELECTRONIC CORRELATION AND 
ENTANGLEMENT 


It was shown by Juhasz and Mazziotti ll35l that the Erobe- 
nius norm of the second cumulant matrix (|pA||^ = Tr[(^A)^]) 
possesses a number of properties that make it a useful mea¬ 
sure. ipAII/r scales linearly with the system size and vanishes 
for single-determinant states. Although it is well suited to 
compare different configurations of the same system such a 
measure is less suited to compare different systems as its up¬ 
per limit is not known. 

Here come informational measures into play, e. g. von Neu¬ 
mann entropy . The entropies based on the 1-RDM have been 
widely studied. Less known are the entropies based on 2- 
RDM. Carlen and Lieb lf36l considered recently the bipartite 
fermionic states and proved several bounds for the entropy 
based on ^D. Eor a Hilbert space of the Fl\ ® FI 2 Hamilto¬ 
nian and corresponding pi 2 density matrix the von Neumann 
entropy can be computed as 


12 - -Trpi2logpi2, 


where the trace is understood in the sense of the tensor product 
Fl\ ®Fl 2 - The 2-RDM of a A-particle fermionic system can be 
considered as a density matrix of a bipartite fermionic state. 
Eor this case the following bounds are known (cf. Carlen and 
Lieb ESI) 


5(pi2)>21nA + G(l), 
2 

2 S1 “ S12 — 


U-Trp2 


> In 


1 - 


( 10 ) 

( 11 ) 


In the following we compare different correlation measures 
and numerically verify the inequalities. 


A. Applications 


All results so far were represented in some abstract one- 
particle basis. In what follows we will focus on molecular sys¬ 
tems with equal number of spin up {No) and spin down (Np) 
electrons. To treat such systems it is convenient to work in 
the closed shell Hartree-Eock (RHE) molecular orbital (MO) 
basis, distinguish spin up and spin down MO states and to 
denote them as i and i, respectively. This implies some addi¬ 
tional symmetries. Obviously, 'Df = and = 0, 

and we have the following blocks for the 2-RDMs; 




kP 




kP 


( 12 ) 


and same holds for the cumulants. Since 2-RDMs in Eq. ( [T^ 
are obtained as scalar products of the CiCj\i//} and c,cj|iA) vec¬ 
tors we will denote these blocks as AA and AB. Other 1- and 
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2-RDMs are obtained by the use of anti-commutation rela¬ 
tions. In order to emphasize spin degrees of freedom we sepa¬ 
rately compute the electronic correlation measures associated 
with A A and AB blocks. 

For the computation of entropies the density matrices must 
be normalized to have trace equal to one. Thus, we introduce 


1 rii 


D': = — 
■' N, 


‘j 


■ '£>',■ 


1 


- 

N,{N, - 1 ) 




1 


^D'l, 

IJ 


(13) 

(14) 


with Tr['D] = = 1, Tr[2D] = 2/; = 1, and 

Carlen and Lieb showed 0^ that the inequalities ( |10|1 l| i 
are saturated for single Slater determinant states. In fact, it 
can be verified that boundary values are achieved for each spin 
channel separately. For single determinant states the cumulant 
matrix ^A vanishes and 2-RDM is given by its unconnected 
components, ^D'j^ - A 'A^. Specifically, = fifki^ij^ki - 

diiSjk), and ^D'jj - fifk^ijSij in the AA and AB channels. The 
first matrix has Na{Na - l)/2 nonzero eigenvalues equal to 2, 
whereas the second one has NaN/} nonzero eigenvalues equal 
to 1. Correspondingly, the entropies in each channel are 5°^ = 
log jNa(Na - 1) = 21ogAa. -H 0(1), and 5^^ = logNaN/j. In 
the next section we will plot the entropies with respect to these 
reference values. 


IV. NUMERICAL ILLUSTRATIONS 


We consider two simple multielectron system, the LiH 
molecule (Figs. |1|2[ ) and the He ring (Figs. and study 
how different correlation measures vary when the system de¬ 
viates from equilibrium geometries. We use our implementa¬ 
tion Il37l[^ of the algorithm by Olsen et al. Il^ based on the 
graphical unitary group approach HOl for performing full Cl 
calculations and subsequent determination of 2-RDMs, cou¬ 
pled cluster and multireference calculations are carried out 



FIG. 1. (Color online) Potential energy surfaces for the stretched LiH 
molecule. 



FIG. 2. (Color online) Full Cl calculation for the stretched LiH 
molecule performed with the 6-311-n-G(2d,2p) basis set, 20 molecu¬ 
lar orbitals are included in a correlated treatment. The electron corre¬ 
lation increases monotonically as the bond is stretched. This is mani¬ 
fested in the correlational energy (a), the entropies (b), the Frobenius 
norm of the second cumulant matrices (c), and the expectation value 
of the z (g) z operator (cf. <[S ig S]> in Eq. jVj or 'Liju SijS in 
Eq. {^) (d). The entropies are given with respect to the single deter¬ 
minant values and S°g. 


with GAUSSIAN 03 program. The systems represent two dif¬ 
ferent scenarios of electronic correlations adopted in quantum 
chemistry. The diatomic molecule is a typical system with 
importance of dynamic correlations. Numerous Slater deter¬ 
minants contribute to the correlation energy, however, one of 
them is dominant. Thus, even for stretched geometries the 
Hartree-Fock (HF) solution represents a valid starting point 
for treating electronic correlations with single-reference meth¬ 
ods HD. This is not so for the He ring for which static cor¬ 
relations are important. Even for slight deviations from the 
equilibrium the wave-function of the system takes a form of a 
sum of several equally significant Slater determinants. It was 
shown by Benard and Paldus ll42l that restricted HF solutions 
are unstable with respect to spin-unrestricted perturbations for 
a wide range of geometries. This invalidates the restricted 
Hartree-Fock approach and also all the single-reference cor¬ 
related methods based on it. At the same time poses an in¬ 
teresting question on whether such spin-instabilities could be 
detected by some correlation measures. 

Let us consider first a simpler case of a single bond break¬ 
ing. It can be easier treated with quantum chemistry meth¬ 
ods as compared, e. g., with double bond breaking in O 2 , C 2 , 
or triple bond breaking in N 2 dimers. However, the physi¬ 
cal mechanisms are the same. LiH possesses only four elec¬ 
trons (in lcr^2cr^ configuration of the ground state) which 

makes it well suited for exact diagonalization studies. For 
these rather small full electron calculations we used the 6- 
311 -i--i-G(2d,2p) basis set yielding a total of 37 basis functions. 
Even with (6s2p)/[4s2p] contracted basis functions for the H 
atom and (12s6p2d)/[5s4p2d] contraction for the Li atom it 
is capable of representing a substantial portion of the corre- 
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FIG. 3. (Color online) Potential energy surfaces for the symmetri¬ 
cally stretched Hg molecule. 



Rh.h(A) 



V3T5 

20 


2 ® 

20 


(d) 



lation energy (Fig. [^. We used the augmented cc-p5z basis 
set for the verification of our results. Stretching the molecule 
(the equilibrium distance is 1.6 A based on a CCSD calcula¬ 
tion) leads to the breaking of a single covalent bond and to an 
increase of the correlation energy (Fig. [^. Our simulations 
indicate similar behavior for the entropies and norms of the 
second cumulant matrices. Notice that the molecule is in the 
singlet ground state and therefore there are only two indepen¬ 
dent spin blocks: - ^A^a, ^Aba - ^Aab- If HCP has the 

electric field vector aligned perpendicular to the molecule’s 
axis (in z-direction) the expectation value of the dipole mo¬ 
ment vanishes ((z) = 0) and we are in the situation analyzed 
above. The correlated part of the survival probability is given 
by the expectation value of the z®z two-particle operator and is 
also computed in AA and AB channels (Fig.|^. Its dependence 
on the interatomic distance is determined by two factors: the 
norm increase of and the matrix elements reduction of the 
dipole operator. The latter dominates the behavior at large 
inter-nuclear separations. 

Let us look now at the manifestly multi-reference system, 
the Hg ring, and investigate how the correlation measures de¬ 
pend on a single geometric parameter, the nearest neighbor 
distance Bh h- Notice, that such symmetric distortion repre¬ 
sents a somewhat artificial situation as the dimerized state pos¬ 
sesses a lower ground state energy. Multireference SCF in the 
subspace of six electrons and six orbitals is capable of recov¬ 
ering a major part of the correlation energy and also correctly 
predicts the asymptotic state of six unpaired electrons residing 
on six independent H atoms. In contrast, RHF method which 
assumes that each molecular orbital is doubly occupied can¬ 
not produce such asymptotic state and fails shortly above the 
equilibrium distance. Even more drastic divergence shows the 
single reference coupled cluster approach. In fact, it does not 
even converge beyond Bh-h “ 3.6 A. On the same Fig. [^re¬ 
sults of full Cl are shown. The potential energy curve runs al¬ 
most parallel to MCSCF, however, has a slightly higher energy 
because not all molecular orbitals (only 30) were included in 
the calculation. Both methods nicely converge towards the 
asymptotic energy of -3 Hartree. We use the augmented cc- 
pqz basis in the (7s4p3d2f)/[5s4p3d2f] contraction. The sur¬ 
vival probability was computed for the case of electric field 


FIG. 4. (Color online) Full Cl calculation for the stretched Hg ring 
molecule performed with aug-cc-pvqz basis set, 30 molecular or¬ 
bitals included in a correlated treatment. The electron correlation 
increases monotonically as the bonds are stretched. This is mani¬ 
fested in the correlational energy (a), the entropies (b), the Frobenius 
norm of the second cumulant matrices (c), and the expectation value 
of the z (g> z operator (cf. ([5 i8> S]) in Eq. jVj or 'ZijuSjjS in 
Eq. 1^) (d). The entropies are given with respect to the single de¬ 
terminant values 5®^ and Thin lines with arrows denote results 
for X = 1 state. Notice, that this state was only computed for larger 
interatomic distances where it becomes indistinguishable from the 
ground state (see the energy plot). Blue and red shading stands for 
the first and second Carlen-Lieb bound (Eq. O). Black and red 
circles denote results for the 3H“ asymptotic state. 


perpendicular to molecular plane, the expectation value of the 
dipole moment is zero ((z) = 0). Correlation energy of this 
system increases steeply as a function of H-H distance. The 
norms of the second cumulant matrix and the entropies be¬ 
have similarly. As in the case of LiH, the expectation value 
of the direct product z ® z is governed by two counteracting 
factors: increasing ^A and reduction of z. At first sight the 
spin instability is not manifested clearly in these correlational 
measures. They show essentially the same dependence (up 
to some scaling) for AA and AB blocks. The difference be¬ 
comes obvious only by comparing with other electronic states 
that have different spin configurations. 

In the asymptotic limits some electronic states of the Hg 
system permit analytical treatment. Because of the vanishing 
electron repulsion, 6H is the asymptotic ground state with four 
possible values of the total spin. At large Bh h the atomic or¬ 
bitals of six H atoms are not perturbed by the interaction with 
other ions or electrons and each electron is in Is state. Thus, 
the density matrices can be conveniently computed in this ba¬ 
sis. Since AL = Np each quantum state of such six electron 
system is given by a linear combination of the following Slater 
determinants \a\a 2 a-iP\fi 2 Pf}, where {a\,a2,a2„P\,P2,P3} is 
a permutation of orbitals 1 to 6 such that ai < a 2 < and 
Pi < P 2 < P 3 , i e- o: and p string are written in the lexico¬ 
graphic order. The spin part of the basis functions is indicated 
by bar for the spin down states. There are in total = 20 
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such determinants and they can be further divided into four 
groups with multiplicities ?c = 1, 3, 5, 7 each containing 5, 9, 
5 and 1 state, respectively. All twenty states are degenerate 
and we consider a particular example of the x - 1 state in 
which all the Slater determinants enter with the same coeffi¬ 
cient. Other spin multiplicity states can be easily constructed 
as explained in Pauncz Il43l . The two-body density matrix is 
built from combinatorics factors and evaluates as follows: 

= —(1 - 5ii), (15a) 

u j‘ 20 

= -^D‘i = -^(1 - 6 ij), (15b) 

ij j‘ 20 ^ 

and vanishes otherwise. Tracing AA or AB blocks leads to the 
same one-particle density matrix 




N„-l 


k=l P k=\ 




IJ- 


(16) 


It is clear that - (cjcJ ) = ^D'j and ^Q'-. = ^D'-. - 0. One 
verifies that the expectation value of the total spin squared 
operator indeed corresponds to the multiplicity ;>i; = 2 ■ 3 -H 1 = 
7: 




(sj=o / V t 





= 12 . 

u 


(17) 


The cumulant matrix is easily evaluated 


^A'/ = -^A'l = (1 - 6a), 

IJ ji 40 ^ 

- —5u, 

IJ 40 40 ^ 

') i i ^ 1 

=-+ — Sii 

ji 40 40 


(18a) 

(18b) 

(18c) 


This yields the norms |pAyi. 4 ||f = and ||^Aab||f = 

The expectation values of z ® z with respect to '^^AA and ^^AB 
are zeros because the dipole transition moments between li 
states vanish due to symmetry as manifested in the selection 
rules. Notice, that for finite distances there is a hybridization 
between these and higher angular momentum states. There¬ 
fore, z ® z can be different from zero. 

The eigenvalues of the matrices (Eqs. ( |T5] l) are is “ 13’ 

36 “ 15 “ ^6 36 “ After the normal¬ 

ization we obtain the pair entropies Saa = Sab - log 15. 
Similarly, the entropies can be computed from the 1-RDM, 
Sa^Sb^ log6. 

Another interesting case is in which the system dissociates 
into three hydrogen anions (H ) and three protons. This is an 
excited 20-fold degenerate, yet bound state of the system with 
the binding energy of just 3 x 0.7542 eV=3 x 0.0277 Haitree 
(measured with respect to the energy of 3 neutral H atoms). 
The state is extremely correlated with RHF not being capable 


of yielding the negative energy (for an overview on the elec¬ 
tronic structure of we refer to A. R. R Rau ll44l . in partic¬ 
ular we quote from this work: ’’The (wave-) function exhibits 
a radial ”in-out” correlation between the electrons such that 
when one electron is ”in” close to the nucleus, the other is kept 
’’out”.”). We can also view as a partial case of the He iso- 
electronic series with the nuclear charge Z being very close to 
the critical value Zcut ~ 0.911 at which ”a quantum phase tran¬ 
sition” from a bound to an unbound two-electron system oc¬ 
curs. A detailed analysis of this system was performed in the 
context of strictly correlated-electrons functional g5l . The 
fact that there are just two electrons in the system allows to 
exactly compute the Kohn-Sham (KS) potential by simple in¬ 
version of the KS equation ll46l . The KS molecular orbitals 
are eigenfunctions of the non-interacting Schrodinger equa¬ 
tion. For the singlet x - I two-electron state there is a one-to- 
one correspondence between the exact density and the Kohn- 

Sham orbitals t//(r) - Fven though the KS poten¬ 

tial can be explicitly constructed, the density functional theory 
does not permit to find the off-diagonal elements of 1-RDM. 
Therefore, we use again the full Cl approach. It is sufficient 
to diagonalize the Hamiltonian for a single H“ atom (for such 
a system the larger aug-cc-pv6z basis set can be used) and to 
construct the for the total system by the observation that 
it is given by a direct sum of the cumulant matrices of each 
subsystem. In order to compute the entropies we reconstruct 
the 2-RMD from the cumulant matrix and using the fact that 
1-RDM is non-zero only when both indices denote states of 
the same H“ anion. 

Numerical results for different correlation measures are 
shown at Fig.|^for the ground state together with asymptotic 
values for the x - 1 state (obtained analytically) and for the 
3H“ state. The latter is based on the full Cl treatment of H“ 
and is only valid for Rh-h For finite distances there is a 

coupling between the subsystems and one needs to diagonal¬ 
ize the full Hamiltonian. This is, however, a formidable task 
since 3H“ is a highly excited state of the system not amenable 
to the Fanczos diagonalization. 

The blue shaded area on the entropy plot shows the differ¬ 
ence between the exact value of the two-particle entropy and 

the first Carlen and Fieb estimate, 25 1 - In ( , ^ , |. The red 

\ i-TiPi / 

shaded area illustrates the second Carlen and Fieb inequal¬ 
ity, 25 1 - ( i_g-si )■ Since at the asymptotic limit the 1-RDM 
and, therefore, 5 1 and Trpj are the same for all twenty de¬ 
generate states, same inequalities hold for all of them, x - 1 
state saturates the inequalities and hence has the largest pos¬ 
sible two-particle entropy among these states. 3 Ha on the 
other hand, has the smallest entropies and Frobenius norms 
out of all states that we considered. For them correlational 
measures lie in between these extremes. In fact, all three of 
them are capable of discriminating between the spin configu¬ 
rations. The survival probability which is closely connected to 
the averaged value of z ® z can even be measured experimen¬ 
tally. In contrast, the 1-RDM and all associated correlational 
measures are ignorant to the spin configuration as the exam¬ 
ple of the asymptotic 6H states shows. Finally, we note the 
correlation between the saturation behavior seen in Fig.Qb,c) 
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when Rh-h is increased beyond « 2.5 A and the flatness at 
these distances of the corresponding potential energy surface, 
depicted in Fig. Even though the electrostatics is less rel¬ 
evant at these distances, the non-local nature imposed by the 
(anti)symmetry on the total wave function persists for the iso¬ 
lated expanding system. This means on the other hand, that 
the asymptotic values of the entropy and |pA||/r delivers short- 
range (atomistic) information. This picture might change in 
a nontrivial way when coupling the electronic and the ionic 
system to external baths, in particular if the bath implies some 
spin-flip scattering. 

V. CONCLUSIONS 

We explored theoretically the possibility of accessing the 
electron correlations using a pump-probe technique with half¬ 
cycle pulses. The non-resonant kick-type excitation is ap¬ 
propriate in sofar as it allows for many-particle excitations 
and can still be amenable to quasi-analytical treatment. We 
demonstrated that to the leading order in the field intensity, 
the probability for a system to remain and/or recur to the ini¬ 
tial state (the survival probability) can be expressed in terms of 


2-RDM and, hence, represents a measure for electronic corre¬ 
lations. Exact diagonalization studies were performed on two 
molecules LiH and Hg - prototypes of systems with impor¬ 
tance from the point view of dynamic and static correlations. 
In the dissociative limit the second system permits for an an¬ 
alytical treatment of some electronic (asymptotic) states for 
which different correlational measures and entanglement en¬ 
tropies were computed. Even though these asymptotic sub¬ 
systems can be considered as noninteracting, the states are 
entangled due to the wave-function anti-symmetry. This is re¬ 
flected in the two-particle entropy, the Erobenius norm of the 
second cumulant matrix ipAUg- and in the expectation value of 
the z®z operator. The latter can be inferred from the proposed 
measurement of the survival probability. Thus corresponding 
experiments can probe the spin configuration of many-body 
systems. 
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